Induced four fold anisotropy and bias in compensated NiFe/FeMn double layers 
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A vector spin model is used to show how frustrations within a multisublattice antiferromagnet such 
as FeMn can lead to four-fold magnetic anisotropies acting on an exchange coupled ferromagnetic 
film. Possibilities for the existence of exchange bias are examined and shown to exist for the case of 
weak chemical disorder at the interface in an otherwise perfect structure. A sensitive dependence 
on interlayer exchange is found for anisotropies acting on the ferromagnet through the exchange 
coupling, and we show that a wide range of anisotropies can appear even for a perfect crystalline 
structure with an ideally flat interface. 



I. INTRODUCTION 

It is interesting to note that some of the technologically 
most important exchange bias systems are also some of 
the most complex and difficult to understand. Antiferro- 
magnetic metal compounds, such as FeMn, can be used to 
pin ferromagnetic layers, making them attractive for ap- 
plication in some devices. The unidirectional and higher 
order magnetic anisotropies that appear in bilayers con- 
taining FeMn seem to vary widely between experiments, 
and there is no model yet capable of explaining the un- 
derlying microsopic mechanisms. 

One particularly intriguing unanswered question is 
how exchange bias can occur in structures with compen- 
sated interfaces. Being compensated, the ferromagnet 
spins in such structures interact equally strongly with 
spins from all antiferromagnctic sublatticcs. In the sim- 
plest approximation, there is no net magnetic moment 
in the antiferromagnet for the ferromagnet to couple to, 
and hence no way for the antiferromagnet to bias a mag- 
netization loop. A closer examination reveals that the 
antiferromagnctic order at the interface is likely to be 
frustrated, and a new configuration resulting in a small 
magnetic moment at the interface should form.^ . However 
it has been shown that this so-called 'spin-flop' coupling 
is not of its own accord able to support exchange bias 
during a magnetization loop measurement^. 

A number of considerations have been discussed that 
may explain the existence of exchange bias in compen- 



sated interface structures. It is highly likely that the in- 
terfaces are not perfectly compensated due to defects in 
structure and chemical composition. These imperfections 
of the interface give rise to small numbers of uncompen- 
sated spins that can result in weak bias shifts. 

A second unanswered question is why the magnitude of 
the bias shift is much smaller than the exchange field cou- 
pling the ferromagnet and antiferromagnet. A number 
of possible explanations have been put forward. These 
include the formation and pinning of a partial domain 
wall near the interface*, the pinning of domains and do- 
main walls in the antiferromagnet^*Sil, and interactions 
between grains in thin antiferromagnctic filmsSi. In each 
explanation there are either restrictions on film thickness 
or reliance on the existence of structural or chemical dis- 
order that make comparison to experiment difficult. 

As most models proposed so far focussed on explaining 
the magnitude of the bias shift, they do not account for 
the drastically increased coercivity observed in exchange 
bias systems. While spin-fiop coupling is not able to shift 
the magnetization loop it does give rise to increased co- 
ercivity in systems with a two sublattice antiferromagnet 
by inducing a two- fold anisotropy in the ferromagnet^. 

In the present paper we demonstrate how frustrations 
within a multisublattice antiferromagnet such as FeMn 
can lead to four-fold magnetic anisotropies acting on the 
ferromagnet. As has been shown previously the inter- 
play between unidirectional and induced higher order 
anisotropies acting on the ferromagnet not only causes 
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an increased coercivity but can also account for some of 
the complex dependencies on applied field angle observed 
in exchange bias systems2iilliiiiiSiiiiiii£. 

The possibility of exchange bias for a multisublat- 
tice antiferromagnet is examined and shown for the case 
of weak chemical disorder at the interface in an oth- 
erwise perfect structure. Most importantly, we iden- 
tify a sensitive dependence for the anisotropics on in- 
terlayer exchange, and show that a wide range of in- 
duced anisotropics can appear even for a perfect crys- 
talline structure with an ideally flat interface. 

The paper is organized as follows. A vector spin model 
for FeMn is introduced in the next section, followed by 
results obtained in the limit of large interlayer exchange 
coupling. Interesting possibilities for multiple configu- 
rations and effective anisotropics appear with small in- 
terlayer coupling, and are discussed in section III. The 
possibility of exchange bias in a perfectly compensated 
system is discussed in section IV. 



II. VECTOR SPIN MODEL AND BULK 
EQUILIBRIUM SPIN STATES 

Magnetic order in metallic antiferromagnets such as 
FcsoMnsois very difficult to predict from first princi- 
ples. One particularly difficult aspect of the problem is 
the importance of contributions that appear as magnetic 
anisotropics that are ultimately determined by spin orbit 
coupling effects inside a crystalline geometry. The prob- 
lem of determining spin configurations associated with 
exchange bias further requires consideration of a large 
number of atoms in both the ferromagnet and antifer- 
romagnet films. These considerations make the problem 
very difficult to approach from a quantum mechanical 
point of view. 

In view of these difficulties, phcnomcnological ap- 
proaches are useful for developing insights into the prob- 
lem. In this work we represent the spin configuration 
using vector spins arranged at lattice sites that represent 
the atomic ordering of a NigiFcig /FcsoMnso two film ex- 
change coupled crystalline structure. Equilibrium spin 
configurations are found as steady state solutions to cou- 
pled sets of classical torque equations for each spin 
located at site i: 



dt 



on 
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(2.1) 



The first term in this equation represents free precession 
of spin Si in its local field calculated from an appropriate 
Hamiltonian, Ti. The second term is dissipative with a 
form chosen to preserve the length of the spin vector. 
The calculation is intended to find equilibria only, and so 
the parameters 7 and a are used only to control stability 
and convergence so that the resulting dynamics does not 
represent specific physical processes. 

The Hamiltonian is chosen to be in the form of a 
Heisenberg spin array with exchange interactions J^.j, a 



Zeeman term for an external applied field H , and four- 
fold anisotropics K^^"^ and -RTf^' . The exact form used 
is 



s 



(2.2) 



The exchange sum is over nearest neighbors at sites i + d, 
and the constants g and hb are the Lande factor and 
Bohr Magneton, respectively. Units are such that Eq. 
12.21 is an energy density. 

All calculations reported here assume fee lattice struc- 
ture for both ferromagnet and antiferromagnet. The 
anisotropy of the antiferromagnet is choosen to mimic the 
3Q phase of an fee antiferromagnet. It is useful to note 
that there are three simple stable equilibrium configura- 
tions, and it is not entirely clear from experiment which is 
favored in FesoMnatjiSiiLiSiiS. Recent calculations sug- 
gest the 3Q phase to be the stable low energy configura- 
tion of bulk Fe5oMn5o^2i2i. The 3Q phase is realised in 
the vector spin model if the anisotropics fulfill both of the 
conditions k\^^ / k\^^ < -1/9 and K^^^ / K^^ < -4/9. 
This phase has spins aligned along cube diagonals as il- 
lustrated in Fig. nia). AH calculations were made using 
periodic boundary conditions in all directions. 

The spin structure calculated using the above model 
for FcsoMnso is strongly affected by surfaces in a thin 
film geometry, due to the broken translational symme- 
try at the interfaces This especially is true when the 
antiferromagnetic film is exchange coupled to a ferro- 
magnet. An example is shown in Fig. ^ (b) where two 
atomic layers of ferromagnetically coupled spins are ex- 
change coupled to an FcsoMuso film. Antiferromagnetic 
anisotropics are k\^'' =—\Jap |/10, and Kf''^ =\Jaf |/10 
where Jaf is the exchange coupling in the antiferromag- 
net. No anisotropy is assumed in the ferromagnet in or- 
der to mimic properties of NigiFcig . The ratio of fer- 
romagnet to antiferromagnet exchange energies is set at 
Jf/\Jaf I = 1-55. These values represent exchange en- 
ergies consistent with the ratio of ordering temperatures 
of NigiFcig and FesoMnso . 

The calculation of a spin configuration is performed in 
analogy to a field cooling procedure. The ferromagnet 
spins are aligned parallel in plane at an angle um made 
with respect to the (001) direction. The antiferromagnet 
spins are initially oriented randomly. The equations of 
motion are integrated numerically until the condition 



y— 



At 



< e 



(2.3) 



is satisfied where e is taken to be on the order of 10"^. 

An interesting feature of the spin configuration shown 
in Fig. n^b) is the near complete antiparallel alignment 
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(a) 



(b) 



FIG. 1: The 3Q phase for bulk FesoMnso are shown in 
(a). An example configuration of spins is shown in (b) for 
FesoMnso exchange coupled to two layers of NigiFeig . The 
dark spheres indicate locations within NigiFeig . Periodic 
boundary conditions are assumed in all directions in (a), and 
in directions parallel to the film planes in (b). 



of spins directly at the interface due to the strong cou- 
pHng between the ferromagnet and the antiferromagnet 
assumed in this case. This leads to a strong modification 
of the entire spin structure of the antiferromagnet to the 
extend that it no longer resembles the 3Q structure. The 
ordering at and near the interface is sensitive to the mag- 
nitude of the interlayer exchange couphng Jf,af , and a 
number of different configurations with similar energies 
are possible. It will be shown below that this has con- 
sequences on the effective anisotropics acting on the fer- 
romagnet, and also on possible mechanisms for exchange 
bias. 



III. INDUCED ANISOTROPY AND 
COERCIVITY 

The magnetic anisotropics in the antiferromagnet are 
communicated to the ferromagnet through the interac- 



tion at their common interface. Cases of strong exchange 
coupling across the interface and weak exchange coupling 
show strikingly different features because of how spins or- 
der in the antiferromagnet. 



A. Strong interlayer coupling 

The strong interlayer exchange responsible for the an- 
tiparallel ordering of spins near the interface in Fig. ^b) 
provides a mechanism for anisotropics to be induced in 
the ferromagnet. If the orientation of the ferromagnet 
is changed through application of an external applied 
field, spins in the antiferromagnet will be rotated through 
anisotropy easy and hard axes. This affects the total en- 
ergy of the system, and in particular can be represented 
by effective fields acting on the ferromagnet. The en- 
ergy per spin is defined for Np ferromagnet spins and 
Naf antiferromagnet spins in a unit cell of the two film 
structure as 



E 



£f.af 



F.AF 



Nf 



Nai 



(3.1) 



where Ef,af is the total energy of the two film system. 
The energy density Ef.af for an Nf = 4 and Naf ~ 12 
system, as depicted in Fig. ^ is shown in Fig. I3a). The 
energy is shown as a function of angle a m and is calcu- 
lated by fixing the orientation of the ferromagnet spins 
and allowing the antiferromagnet spins to relax using the 
field cooling procedure described earlier. The field cool- 
ing process was repeated for each orientation value of 
Q!M • The results of this procedure were compared to 
results from an alternate method by which the antifer- 
romagnet was field cooled along um — and the energy 
calculated for each new um without additional field cool- 
ing. The two different calculation procedures produced 
identical results. 

It is clear from Fig. E^a) that interlayer coupling be- 
tween the ferromagnet and antiferromagnet results in a 
total energy for the system with four-fold symmetry for 
the in-plane orientation of the ferromagnet. The solid line 
in Fig. 121 is a fit using 



£f,af 



Kl^^ sin^(aM - 04) cos^(q;a/ —04) (3.2) 



where describes the orientation of the anisotropy easy 
and hard axes with respect to the [100] crystallographic 
direction. 

Magnetization loops calculated for this set of K^^^^ , 

Kl , and Jf.af parameters show properties consistent 
with a simple four fold anisotropy of the form in Eq. 
13.21 An example is shown in Fig. Olb) for a field ap- 
plied along the [110] and [100] directions. The magne- 
tization is given in units of Bohr Magneton and is de- 
fined as Us = —gfJ'B X^r^ I Si I- The easy direction is for 
a4 = 7r/4. There is coercivity for the field along the easy 
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FIG. 2: Total energy per spin is shown in (a) as a function 
of ferromagnet orientation. The angle cvA/is measured with 
respect to the [100] axis. The effective anisotropy displayed by 
the ferromagnet is four fold, without a bias shift. An example 
magnetization loop is shown in (b) for the applied field along 
the [110] (filled circles) and [100] (solid line) directions. 

direction as would be expected due to the existence of 
stable and metastable states for the ferromagnet parallel 
and antiparallel to the applied field. Coercivity along the 
hard direction exists because of the zero field remanent 
magnetization aligned along an easy direction. 



B. Moderate and weak interlayer coupling 

In the case of large Jf.af the magnetic behaviour of 
the ferromagnet is dominated by the intra- and interlayer 
exchange coupling and is well described by a simple four 
fold anisotropy. When Jf,af is not large relative to Jaf , 
new possibilities for metastable ordering within the anti- 
ferromagnet appear. Examination of results from numer- 
ical solutions to Eq. 12. II reveal multiple equilibrium con- 
figurations with comparable but different energies. An 
example of this behavior is shown in Fig. |31 where the 
the total energy ef.af^s shown as a function of angle 
for an interlayer exchange coupling Jf.af = —0-5\Jaf \- 
The parameters and geometry are otherwise as used for 
the example shown in Fig. |21 These results were gener- 
ated using the field cooling procedure at each angle. At 
each angle, thirty different random initial configurations 
for the antiferromagnet spins were used, resulting in a 
spread of energies as shown in Fig. Ota). The range of 



energies at each angle represents a sampling of different 
possible spin configurations. 

The open symbols connected by the thick solid line in 
Fig. ISJa) are the lowest energies found during the cool- 
ing process. In Fig. OJb) the lowest energies are plot- 
ted separately to more clearly show the four-fold symme- 
try of the ground state energies. Note that the energy 
eF,AF displays very sharp maxima and is only approxi- 
mately consistent with the K^-^^ as given in Eq. 13.21 



(a) 
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FIG. 3: Total energy as a function of angle with Jf,af ~ 
—0.5\Jaf |. Other parameters and structure are the same as 
used for Fig. |5| In (a), energies from thirty different ran- 
dom initial configurations calculated at each angle are shown. 
The spread in energies at each angle is due to the existence 
of several metastable spin configurations within the antifer- 
romagnet. The lowest energies found are identified by open 
symbols in (a) and plotted together in (b). 

The degree to which the lowest energy configurations 
result in an anisotropy approximating that described by 
Eq. I3.2l is strongly dependent on the strength of the inter- 
layer coupling Jf,af ■ For some values of Jf,af , even the 
sign of the effective anisotropy can change, representing a 
change in the orientation of the easy and hard axes. The 
magnitude and sign of K^'^'^ saturates to a fixed value 
for Ji?_AF larger than 1.25\Jaf \ - At these large values, 
the antiferromagnet spins at the interface are aligned 
coUinear with respect to the ferromagnet spins, and ro- 
tate rigidly with the ferromagnet as described in the pre- 
vious chapter. 

A plot of the effective anisotropy determined by fit- 
ting the ground state energy configurations to Eq. 13.21 
is given in Fig. 0] as a function of interlayer coupling 
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Jf,af ■ The goodness of fit is quantified by the measure 
where = 1 is very good and denoted by Ughtly 
shaded stripes, and very bad bits are denoted by dark 
gray stripes. 



(a) 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



FIG. 4: Effective anisotropy describing ground state energy 
configurations as a function of interlayer coupfing Jf,af ■ The 
degree to which the anisotropy is a simple four-fold type of the 
form given in Eq. I3.2l is indicated by grayscale shading. White 
means that the description is very good, and black means that 
the description is very poor. Note that the direction of the 
easy and hard axes can change depending on the magnitude 
of Jf,af ■ 

It is interesting to examine separately the ferromagnet 
and antiferromagnet contributions to Ef.af ■ In particu- 
lar, the ferromagnet component of the energy can have 
a different dependence on orientation angle a m than the 
total system energy. A way of thinking of this is to con- 
sider the effective field he// acting at a layer of spins 
within the structure and calculate the associated energy 
—Si • he//. The field varies in magnitude and direction 
throughout the structure, and can have a dependence on 
um that also varies from layer to layer. A way of charac- 
terizing this difference in a meaningful way is to identify 
extrema in the energy of the ferromagnet calculated as a 
function of aM , and compare this to the extrema deter- 
mined from the total energy e f,af ■ 

Results of this characterization are given in Fig. |31 as 
the angles aM where minima (a) and maxima (b) occur as 
functions of Jf,af ■ The angles amin,F represent orien- 
tations of the ferromagnet spins where the energy of the 
ferromagnet spins has minima, and the angles am,in,F+AF 
represent orientations of the ferromagnet spins where the 
energy of the entire system has minima. Likewise, the 
angles amax,F represent orientations of the ferromagnet 
spins where the energy of the ferromagnet spins has max- 
ima, and the angles amax,F+AF represent orientations of 
the ferromagnet spins where the energy of the entire sys- 
tem has maxima. 

An interpretation of the minima is to assign different 
easy and hard axes to the ferromagnet and system. In 
this view, the alignment of the ferromagnet to an easy 
axis associated with the system rather than the easy axis 
defined by the local effective field acting on the ferromag- 
net is a clear indication of how order in the antiferromag- 




- O O 



0.5 0.0 




FIG. 5: Orientations for energy extrema as functions of 
Jf,af ■ In (a), angles ami'n,F represent orientations of the 
ferromagnet spins where the energy of the ferromagnet spins 
has minima, and the angles amin,F+AF represent orientations 
of the ferromagnet spins where the energy of the entire sys- 
tem has minima. In (b) the angles amax.F and amax,F+AF 
are shown for the corresponding energy maxima. 



net is involved in determining the magnetic anisotropics 
observable through the ferromagnet. As an extreme ex- 
ample, the large Jf,af hmit shown in Fig. |5l has the 
ferromagnet 'easy' axis parallel to the system hard axis. 
This occurs again at smaller values of Jf,af ■ 




FIG. 6: Coercive fields for applied fields along the [110] and 
[100] directions for different interlayer couplings Jf,af ■ The 
closed symbols are for the applied field directed along the 
[110] and the open symbols are for the applied field directed 
along [100]. 

Coercivities calculated for magnetization loop Simula- 
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tions also show curious behavior for some values of in- 
terlayer coupling. Coercive fields for applied fields along 
the [110] and [100] directions are shown in Fig. El The 
most curious feature is the different coercivities produced 
with positive and negative Jf,af ■ The reason for this is 
that the antiferromagnetic film used for the calculation 
is very thin. Large values of \Jf.af \ fuhy align the spins 
colinear with the ferromagnet and contribute to the net 
magnetic moment. When Jf,af > 0, the contribution 
increases the Zeeman energy of the system in an applied 
field, and when Jf,af < 0, the contribution decreases 
the Zeeman energy in an applied field. The easy axis 
coercivities are therefore larger with negative Jf.af than 
with positive Jf,af ■ 

IV. MECHANISMS FOR EXCHANGE BIAS 

The anisotropics induced on the ferromagnet through 
interlayer exchange coupling with the ferromagnet were 
not found to contain any unidirectional contribution in 
any of the compensated interface examples studied. This 
is consistent with previous calculations demonstrating 
the inability of spin flop coupling at compensated in- 
terfaces to support exchange bias shifts using physically 
reasonable assumptions for anisotropy field a^'^^ . 



A. Bias with interface defect 

It has been noted that small regions of uncompensated 
spins at the interface can be sufficient to create exchange 
^jg^g3j7j23j24 'pj-^jg possible mechanism for exchange bias 
is explored for the NigiFcig /Fe5oMn5o model discussed 
here. The unit cell of the structure with defect is based 
on eight antiferromagnet atomic layers with eight spins in 
each layer exchange coupled to a ferromagnet film con- 
sisting of two atomic layers, also with eight spins per 
layer. Periodic boundary conditions in the planes paral- 
lel to the interface as used. The defect is represented by 
replacing the exchange couplings and local anisotropics of 
one spin on the ferromagnet side of the interface with val- 
ues appropriate to the antiferromagnet. In this way the 
unit cell of the structure has 15 ferromagnetically cou- 
pled spins, and 65 antiferromagnetically coupled spins, 
with a small uncompensated region at the interface. 

Results from calculations for coupling between fer- 
romagnet and antiferromagnet spins with Jf.af — 
— 0.5] J^F 1 are shown in Fig. [7| A small negative bias of 
magnitude 0.005] Jaf Saf\/ (giJ-s) appears. It is interest- 
ing to note that a simple estimate of the shift would be to 
average the uncompensated coupling energy 2 Jf.af over 
the number of ferromagnet spins, giving a bias field of 
0.07\ Jaf SAF\/{gfJ'B)- The factor of fourteen discrep- 
ancy is because the spin order in the antiferromagnet 
changes during magnetization in such a way as to reduce 
the magnitude of the field necessary to align the ferro- 
magnet. This is analogous to the reduction in bias field 





1 ' 























-0.02 -0.01 0.00 0.01 0.02 



FIG. 7: Biased magnetization loop with one defect spin. The 
parameters used are the same as those used for the unbiased 
loop shown in Fig. |Hb). The difference is that here the 
unit cell has Nf = 15, Naf ~ 65, and one spin 'defect' 
located at the interface providing a small net uncompensated 
antiferromagnetic moment. Jf.af is —0.5\Jaf \- 

described by Mauri, et al^ for exchange bias with com- 
pletely uncompensated interfaces. 

V. SUMMARY 

A molecular field model of magnetic order in multi- 
sublattice magnets, such as FcsoMnso , exchange coupled 
to a ferromagnet has been examined. It has been shown 
that ordering of spins near the surfaces and interfaces 
of the antiferromagnet are strongly affected by exchange 
coupling, but that the intrinsic four-fold anisotropy of 
the antiferromagnet is still induced into the ferromagnet. 
Through exchange coupling, the ferromagnet spins expe- 
rience a four-fold anisotropy with a magnitude sensitive 
to the strength of the interlayer coupling. An important 
point is that the nature of the anisotropy is also depen- 
dent on the strength of the coupling, and is only strictly 
of a simple four fold form in the case of strong interlayer 
coupling. 

Weak and moderate values of interlayer coupling 
lead to an additional interesting effect on the induced 
anisotropy. Rather than a single well defined order in 
the antiferromagnet, a number of metastable configura- 
tions appear. The lowest energy configurations lead to 
an induced anisotropy with four-fold symmetry. The de- 
gree to which the four-fold anisotropy associated with the 
lowest energy configurations is described by a simple si- 
nusoid depends on the strength of the interlayer exchange 
coupling and is only approximate at best for some values 
of the exchange. 

Exchange bias shifts for perfectly compensated inter- 
faces assuming the low energy spin configuration was not 
found. Instead, an exchange bias could be created by 
allowing the interface to have some small degree of un- 
compensation through introduction of a defect in the in- 
terface structure. These findings are consistent with pre- 
vious work showing that a spin-flop configuration at a 
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compensated interface is incapable of producing a shifted 
magnetization loop, whereas mixed interfaces with some 
amount of uncompensation present can shift the magne- 
tization loop. 
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